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We consider a five-dimensional Minimal Supersymmetric Standard Model compactified on a S 1 /Z2 
orbifold, and study the evolution of neutrino masses, mixing angles and phases for different values 
of tan P and different radii of compactification. We consider the usual four dimensional Minimal 
Supersymmetric Standard Model limit plus two extra-dimensional scenarios: where all matter su- 
perfields can propagate in the bulk, and where they are constrained to the brane. We discuss in both 
cases the evolution of the mass spectrum, the implications for the mixing angles and the phases. We 
^ ' find that a large variation for the Dirac phase is possible, which makes models predicting maximal 

, leptonic CP violation especially appealing. 

\0 ■ I- INTRODUCTION 

(N _ 

In the last decade, our perspectives in the search for physics Beyond the Standard Model (BSM) have seen the 
development of theories with extra (compact) dimensions. TeV-scale extra-dimensional models have allowed us to 
Oh build effective models which can be tested at the present and next generation of colliders, where these wide varieties of 
extra-dimensional models have been proposed to solve, or at least understand from a geometrical perspective, different 
. theoretical problems arising in the Standard Model (SM) and in its four-dimensional space-time extensions. 

Another ingredient which can be introduced, supersymmetry, and in particular its minimal low energy construction, 
the Minimal Supersymmetric Standard Model (MSSM), has long been held to play a central role in BSM physics; 
even if at present there is no evidence at colliders for the supersymmetric partners of the SM particles. From a 
\ theoretical point of view, supersymmetry plays a key role in resolving many problems in the SM and BSM physics; 
00 ' from gauge coupling unification, to the hierarchy problem, to the building of realistic grand unification models. The 
combination of supersymmetry and the physics of extra-dimensions has the added bonus of stabilising the extra- 
dimensional theory from quantum fluctuations, as well as the extra-dimensions potentially providing a mechanism 
for supersymmetry breaking. Indeed, four-dimensional supersymmetric models typically lack a simple mechanism for 
supersymmetry breaking, which the extra-dimensions may offer. We will discuss in the following a five-dimensional 
Af = 1 supersymmetric model compactified on the S 1 /Z2 orbifold (5D MSSM) as a simple testing ground for the 
effects of the extra-dimension on the neutrino sector. A similar study, concerning the quark Yukawa couplings and the 
CKM matrix observables was performed in Ref. [l[ , which contains also a more detailed presentation of the general 
features of the 5D MSSM model. 

Neutrinos are generally taken to be massless in the SM, and the experimental evidence for nonzero neutrino masses 
implicit in the neutrino oscillations measurements, gives an important indication for physics BSM. In fact, neutrino 
masses are many orders of magnitude smaller than those of quarks and charged leptons. However, contrary to the 
small mixings in the quark sector, two of the lepton mixing angles are identified as being rather large, close to maximal. 
For an overview of the present knowledge of neutrino masses and mixings see Ref. [Jul an d references therein. The 
most recent experimental evidence of this fact is the measurement of the 813 mixing parameter by the Daya Bay 
and RENO experiments [1, The implication of three sizable mixing angles are huge, and will surely boost the 
number of investigations of the neutrino mixings and phases in the near future; in particular, this result suggests the 
possibility of measuring leptonic CP violation. The neutrino sector seems, therefore, to continue to play a special role 
in understanding BSM physics. 

Recall also that quark and lepton masses and mixing angles are free parameters in minimal extensions to the SM. 
The running of both quark and neutrino masses and mixings has been investigated extensively in the SM and various 
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BSM extensions Ul|6Hl5|. Also, as neutrino mixing angles show a pattern that is completely different from that of 
quark mixings, the relative wealth of the latest experimental data has motivated efforts on the theoretical side to 
understand the possible patterns of neutrino masses and mixings, and therefore to expose the underlying fundamental 
symmetries behind them. As such, understanding the evolution of these neutrino sector parameters will be critical as 
higher energy experiments probe this sector. With this in mind, starting from the 5D MSSM, this paper is organized 
as follows: In section [TT] we introduce the dimension five operator (in four space-time dimensions) for the neutrino 
masses in the 5D MSSM. Then in section 11111 we further develop the renormalization group equations for the neutrino 
mass matrix. In section IIVI we shall explore and discuss the evolution properties and behaviours for neutrino masses 
and mixing angles, and section fVl is devoted to our conclusions. 



II. 5D MSSM AND NEUTRINO MASSES 



In order to study qualitative features in a model independent way, an attractive and simple possibility is to use 
a low-energy effective theory formulation. This means that additional particles and symmetries present at higher 
energies are organized within a systematic expansion. We shall consider the heavy states arising from BSM physics to 
decouple at low energies. In that case, the degree of suppression of an operator in the low energy effective Lagrangian is 
characterized by its mass dimension. This general framework is valid also when discussing extra-dimensional theories. 
In the 5D MSSM, the Higgs superfields and gauge superfields always propagate into the fifth dimension. However, 
different possibilities of localisation for the matter superfields can be studied. We shall consider the two limiting cases 
of superfields with SM matter fields all in the bulk or all superfields containing SM matter fields restricted to the 
brane. When all fields propagate in the bulk, the action for the matter fields <&i is [T^ ]: 

Smatter = J d* zd V + + ^8^8(6) - ^8^6(6) 

+5(2$^ - 2^v^ + $tx$Ad) + $ixms(0))} ■ (i) 

Similarly, when all superfields containing SM fermions are restricted to the brane, the part of the action involving only 
gauge and Higgs fields is not modified, whereas the action for the superfields containing the SM fermions becomes: 

Smatter = J d s zdyS{y) + 2g$ t V<S> z } . (2) 

Neutrinos, being SM matter fields, will therefore be considered either localised or propagating in the bulk (in these 
two different scenarios). However, as the neutrino masses are roughly six orders of magnitude smaller than the other 
light SM fermions, and if neutrinos are Majorana particles, such a small mass could be understood if there is new 
physics beyond the electroweak scale. The lowest order operator, which generates Majorana neutrino masses after 
electroweak symmetry breaking (EWSB), is the lepton- number violating Weinberg operator [l7|. This lowest order 
operator (appearing with dimension d — 5 in four space-time dimensions) can be written as: 

-^(L c Je^^)(4e s ^ 7 ) + h.c, (3) 

where L and <j> are the lepton and the Higgs doublet fields. M is the typical heavy energy scale for the range of 
validity of the low-energy effective theory. An operator of this type can be generated, for instance, by the usual see- 
saw mechanism. In which case the heavy scale M can be identified with the mass of the heavy right-handed neutrino. 
After EWSB the Higgs acquires a vacuum expectation value (vev) and the operator in Eq. ([3]) gives a Majorana mass 
term for the neutrinos. In the context of the MSSM it can be written in the form: 

-^.(L^H^(L^H-), (4) 

where L and H u are the lepton and up-type Higgs doublet chiral superfields respectively. This operator is crucial 
for the study of neutrino masses and mixings, where renormalisation group equations for this effective operator have 
been derived in the context of the four-dimensional SM and MSSM pjj, LLJJ • 

For J\f = 1, D = 4 supersymmetry the superfield formalism is well established: superfields describe quantum fields 
and their superpartners as well as auxiliary fields as a single object, thus simplifying considerably the notations and 
the calculations. A similar formulation for a 5D vector superfield has been developed in Ref. [Hj], and the superfield 
formulation for vector and matter supermultiplets can be found in Refs. [l9|, |2Q |. A detailed description of these 
techniques in the covariant formalism, together with the corresponding Feynman rules, is described in Ref. [21] . 
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In Eq. (|4} the neutrinos acquire Majorana masses through the dimension five operator. An extension to the 
compactified 5D MSSM was considered in Ref. [l6j], and we shall use a similar formalism here. That is, for the energy 
range below the mass scale of the heavy right-handed neutrino, we assume that the effective theory is just the effective 
5D MSSM with S 1 /Z 2 compactification, and the energy dependence of the effective neutrino mass matrix below that 
scale is then described by its renormalization group equation. Note that in our model the Yukawa couplings in the 
bulk are forbidden by the 5D TV = 1 supersymmetry. However, they can be introduced on the branes, which are 
4D subspaces with reduced supersymmetry. We will write the following interaction terms, called brane interactions, 
containing Yukawa-type couplings: 

Strane = J d 8 zdyS(y) j |^ A, „v <i vl >.,<!•/. - ^LiHvLjH^J 5(6) + /i.e. | . (5) 

The last term corresponds to the effective neutrino mass operator, with dimensional coupling fey whose evolution and 
(that of the Majorana mass term for neutrinos) we shall now calculate. 

III. k COUPLING EVOLUTION 

In the case that the neutrino superfields propagate in the bulk, it is straight forward for us to write down the 
following terms related to the neutrino masses, where we include terms related to both Majorana and Dirac masses: 

feij . n TT n , n . . 1 1 fei 



AM 



L\BlL)Hl + . 2V2Z»#°£X + h.c. + ± (l*0 L « + L^L™) . (6) 



In which case the first two terms are from the F terms of the dimension five operator in Eq. ([5]) along with its Kaluza- 
Klein expansion in the bulk space, and the last term is from the F term related to the gauge covariant derivative in 
Eq. (jTJ) . i? is the compactification radius of the 5D MSSM. 

In the present case we consider the effective neutrino mass operator with dimensional coupling fey ; after spontaneous 
symmetry breaking, the Majorana neutrino masses can be written as m v = kv 2 sin 2 /? (v being the vev of the Higgs 
field and tan/?, the ratio of the vevs of our two Higgs doublets) and k — fey / '(2MitR) for bulk propagating, and 
k = k/(2M) for brane localised matter superfield scenarios respectively. For simplicity we have the mass matrix for 
one generation of neutrinos as follows: 



( m v 2V2m v ^ 
2V2m u 



v ° s V 



(7) 



with the basis being v T = (^l, )■ By using the inverse iteration method, and up to first order corrections, 

we observe that the three active neutrino flavour eigenstates v a are related to their mass eigenstates i>i through the 
transformation v a « Upmns(1 + 0(i? 2 r7i^))j/;, where we have chosen to work in the basis where the charged lepton 
Yukawa coupling matrix is diagonal, and the neutrino mixing matrix Upmns, up to a diagonal phase matrix, stems 
from the diagonalisation of the Majorana neutrino mass matrix m„. That is 

UpMNS m v u PMNS = diag(m 1 ,m 2 ,m 3 ) , (8) 

with rrii being the minimal extension SM neutrino masses. For the large hierarchy between the compactification 
scale 1/R (which is usually taken at the TeV level for LHC physics) and the neutrino mass (eV level), we conclude 
Va, = UpMNS^i, up to a correction of order 10~ 24 . 

As in our previous works, the evolution of the Yukawa couplings were derived using standard techniques (see for 
example Appendix IBI and Ref. [16|). The evolution of k can then be determined for the two scenarios we wish to 
consider along the same lines as we have done earlier [lfij]. As such, the beta function of k at one loop for the MSSM 
is: 



(167r 2 )/3 fc = ak+([Y c T Y e *]k + k[Y^Y e ^jC 

= (&Tr{Y^Y u ) - - 6. 9 ^ k + {[Y e T Y e *]k + k[YjY e ]j 1 



(9) 
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where C = 1 and a = 6Tr(Y^Y u ) - §</? - 6g 2 = 6(y t 2 + y 2 c + y 2 u ) - f g 2 - 6 5 2 . 

We will use S(fi) = /j,R, t = ln(^-) and 5(i) = e l MzR, where at one loop for the case in which matter fields 
propagate in the bulk: 

(167r 2 )/3 fc - ak + (\Y e T Y e *]k + k[YjY e ty C 

= (W^TMyJiy - (|g 2 + 6g 2 )^ k + (\Y e T Y e *]k + k[YjY e f) tt^R 2 , (10) 

where in this case C = irS{t) 2 and a = 6?:S(t) 2 Tr{Y^Y u ) - (| 5 2 + 6g%)S(t). 

The beta function of k at one loop for the brane localised matter fields case is then 

(167r 2 )/5jt = ak + ([Y e T Y e *]k + k[YjY e ^j C 

= ^6Tr(rJy„)-(^ 5 2 + 9 3 2 ) M i?)A : +([Y; T y;*]fc + fc[y; t y;])2 Ai i?, (11) 

where now we have C = 2S(t) and a = QTr(Y^Y u ) - (§£ 2 + 9g$)S(t). 

From Eq. §E§ and by inserting this into the beta function for k, one obtains the evolution for the neutrino mixing 
matrix due to the Kaluza-Klein modes. Furthermore, in order to derive the renormalization group running behaviour 
of the neutrino mixing parameters, we employ the standard parameterization, in which Upmns is parameterized by 
three mixing angles and three CP-violating phases, see appendix [3] for details. These equations are combined to 
arrive at the runnings for the leptonic mixing angles, as given in appendix [Cj and where we have also included the 
phases in appendix IDl 

Note that as a general remark, it is typical that when all fields are in the bulk that extra-dimensional effects are 
more important in a number of observables and, at the same time, the typical cut-off of the effective theory is lower 
with respect to the brane localised field model (as perturbative unitarity is violated at a lower scale and the strong 
running of many observables may drive the model outside the domain of perturbation theory for which the evolution 
equations can be applied). We will examine in detail the behaviour of the neutrino sector observables with this 
constraint in mind in the next section. 



IV. NUMERICAL RESULTS 



In this section we will apply the beta functions to the study of the evolution of the different couplings. Unless 
explicitly mentioned we refer to the normal hierarchy for the neutrino masses m 3 > rn 2 > mi. As there is some 
freedom in the choice of the mass values for the initial conditions, within the allowed range of Am 2 ol and Am 2 (m , we 
have explicitly used in the case of a normal hierarchy the values for initial low energy conditions indicated at the end 
of Appendix lAl 

The discussion of the fixed point structure for the mixing angles in the infrared is the same as the MSSM one 
[HI [HI i where the evolution equations for the neutrino masses and mixings are formally similar to the MSSM case 
since the beta functions have a similar structure. As expected tan /3 plays an important role as all the mixing angles 
and phases depend on y T . However, the new degrees of freedom (the extra-dimensional fields giving rise to Kaluza- 
Klein excitations of the zero modes) become important at energies corresponding to their masses. In the following 
we study the evolution of the relevant parameters, such as Am 2 ol , Am 2 tm and the angles and phases, as a function 
of the energy scale and of tan ft. Only some selected plots will be shown and we will comment on the other similar 
cases not explicitly shown. 

In general, in the brane case, the evolution has the same form for the three masses mi, m 2 , m^. This leads to a 
reduction of up to a factor of two for the masses at t = 6 (for a large radius, R^ 1 = 1 TeV) with respect to the MSSM 
values at low energies (smaller radii give a weaker effect as the Kaluza-Klein excitations contribute to the evolution 
equations at higher energies) . This prediction is extremely stable and can be explained as the evolution of the masses 
is governed by the equation 

™* = HT^2 mi ( a + Ci Vt) j ( 12 ) 

lD7T z 

where the coefficients Cj induce a non-universal behaviour and the parameter a contains the up Yukawa couplings 
and the gauge coupling terms (detailed in Appendix [C]). In contrast to the MSSM, the evolution in the brane case is 
completely dominated by the universal part. The essential point is that in the MSSM the positive contribution to a, 
approximately 6y 2 , is of the same order as the negative contribution from the gauge part. In our case the gauge part 
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has a large pre-factor S(t) = e t MzR with respect to the MSSM which makes it completely dominant compared to 
any other contribution. As energy increases, we can write: 



1 



16tt 2 



< o 



(13) 



From this approximation we immediately see that all masses decrease with increasing energy and eventually tend to 
zero if the evolution equations can be trusted up to a high energy. 

In the bulk case the evolution has again the same form for the three masses mi, m,2, m^, but the behaviour is the 
opposite as the masses increase at high energy because all matter fields propagate in the bulk and contribute to the 
evolution. In detail this can be seen by the fact that even if the gauge part gets a large pre-factor S(t), the Yukawa 
part gets in this case a pre-factor S(t) 2 , which changes the sign of the derivative with respect to the previous case: 



1 



16-7T 



;rrii 



S(t) 



> 



(14) 



From this approximation we see that all masses increase with increasing energy scale. 

The situation is more involved when analyzing the mass squared differences. We plot in Figs. [T] and [5] the evolution 
of Amj ol and /S.m 2 atm both for the matter fields on the brane and for all fields in the bulk. In the brane case different 
behaviours as a function of the energy scale are possible as a relatively large interval in energy range is allowed for the 
effective theory. As explicitly illustrated in Fig. [TJ the relevant radiative corrections controlled by the gauge fields in 
Eq. (fT3"|) become dominant as energy goes up, which tends to reduce mass splitting, and an approximately degenerate 
neutrino masses spectrum at the high energy scale mi ss m-2 « 7713 becomes favourable. This is in contrast with 
the MSSM, where the neutrino mass splitting becomes large at an ultraviolet cut-off. Therefore, it is very appealing 
that the neutrino mass splitting at low energy could be attributed to radiative corrections resulting from a degenerate 
pattern at a high energy scale. In Fig. [21 the bulk case tends to a non-perturbative regime, where the unitarity bounds 
of the effective theory are reached much faster and only a much shorter running can be followed using the effective 
theory. As seen in Eq. (fT4")l , the quadratic terms related to S(t) dominate during the fast evolution. As such, the 
neutrino mass splitting becomes even larger at a high energy scale. 




FIG. 1: Evolution of Am^j (left panel) and Am^tm (right panel) as a function of the parameter t — ln(/i/Mz) with matter- 
fields constrained to the brane. In this plot we have taken tan/3 — 30. The black line is the MSSM evolution, the red (small 
dashes) is for R^ 1 ~ 1 TeV, the blue (dash-dotted) R' 1 ~ 4 TeV, the green (large dashes) R" 1 ~ 15 TeV. 

Concerning the evolution of the mixing angles, as can be seen in Figs. [3] [H the largest effect is for O12, with changes 
of more than 70% possible for the brane localised matter field scenario. As observed, due to the large quadratic term 
of S(t) in the beta function, the #12 has a rapid and steep variation in the bulk case. However, for the brane case, it 
has a relatively longer evolution track with the O12 then being pulled further down until the termination point (where 
the effective theory becomes invalid). In contrast, the running of #13 and 623 is much milder. As demonstrated in 
Figs. [3] and |31 changes in the values of #13 vary only a couple of degrees. For a larger value of tan/3 we have a 
relatively large Yukawa coupling to r, which leads to a large magnitude for its beta function, resulting in a relatively 
large variation during the evolution. However, a running to #13 = cannot be observed in any situation. From the 
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FIG. 2: Evolution of Am 2 ao i (left panel) and Am^tm (right panel) as a function of the parameter t — ln(p/Mz) with matter 
fields in the bulk. In this plot we have taken tan/? = 30. The black line is the MSSM evolution, the red (small dashes) is for 
R~ x ~ 1 TeV, the blue (dash-dotted) R~ 1 ~ 4 TeV, and the green (large dashes) R~ 1 ~ 15 TeV. The evolution is towards large 
values at high energies, the divergence signaling the breakdown of the effective theory. 

evolution behaviour of #13, one can see that the renormalization group running effects or finite quantum corrections 
are almost impossible to generate #13 = at a high energy scale, even though the power law enhanced evolution is 
considered during the running. Therefore, for the tri-bimaximal mixing pattern [22], in the current context with no 
other extreme conditions being taken into account, a slightly changed #13 could not be accommodated during the 
whole range of the energy scale. Similar trajectories are also observed for #23. 
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FIG. 3: Evolution of 613 as a function of the parameter t = ln(fi/Mz) with matter fields in the bulk. In this plot we have taken 
tan/3 = 5 ( left panel) and tan/3 = 30 (right panel) respectively. The black line is the MSSM evolution, the red (small dashes) 
is for R" 1 ~ 1 TeV, the blue (dash-dotted) R~ x ~ 4 TeV, and the green (large dashes) R~ x ~ 15 TeV. 



Noting that the Dirac phase 5 determines the strength of CP violation in neutrino oscillations. The runnings we 
include follow the general features presented in Figs(7] and [51 with large increases possible once the first Kaluza- 
Klein threshold is crossed. From these studies we have seen that the variation is bigger for high tan/3 with large 
changes appearing in the brane case when approaching the high energy scale. The recent results from the Daya 
bay and RENO reactor experiments have established a non zero values of #13. Therefore, the leptonic CP violation 
characterized by the Jarlskog invariant J ~ sin #12 cos #12 sin #23 cos 623 sin#i3COS 2 #i3 sin 5 becomes promising to be 
measured in the future long baseline neutrino oscillation experiments. As plotted, we can observe a relatively large 
evolution for the Dirac phase, even the maximum CP violation case <5 = f could be achieved for relatively small input 
values. For leptogenesis related to the matter-antimatter asymmetry, we should note that the parameters entering 
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FIG. 4: Evolution of 9i3 as a function of the parameter t = ln(/i/Mz) iviift matter fields constrained to the brane. In this plot 
we have taken tan/3 = 5 ( left panel) and tan/? = 30 (right panel) respectively. The black line is the MSSM evolution, the red 
one (small dashes) is for _R _1 ~ 1 TeV, the blue (dash-dotted) R' 1 ~ 4 TeV, and the green (large dashes) R' 1 ~ 15 TeV. 




FIG. 5: Evolution of Q\i in the bulk (left panel) and on the brane (right panel) as a function of the parameter t = hi(fj,/Mz)- 



In this plot we have taken tan/3 = 30. The black line is the MSSM evolution, the red one (small dashes) is for R 
the blue (dash-dotted) J? -1 ~ 4 TeV, and the green (large dashes) R' 1 ~ 15 TeV. 



1 TeV, 



the leptogenesis mechanism cannot be completely expressed in terms of low-energy neutrino mass parameters. Note 
that in some specific models the parameters of the PMNS matrix (which contains CP asymmetry effects) can be used 
[1,[23[. Here, the CP-violating effects induced by the renormalization group corrections could lead to values of the CP 
asymmetries large enough for a successful leptogenesis, and the models predicting maximum leptonic CP violation, or 
where the CP-violating phase S is not strongly suppressed, become especially appealing. Sp ecific models with large 
extra dimensions in which leptogenesis is relevant at low scale can also be found in Ref. |24| . 

As seen in Appendix [C] the running of the mixing angles are entangled with the CP-violating phases, and a more 
general discussion is beyond the scope of this paper. The phases <pi an d 4>i do not affect directly the running of the 
masses, while the phase S has a direct effect on the size of dm/dt, although its importance is somewhat reduced by 
the magnitude of 6*13. For further discussions of the correlation between these phases and mixing angles, refer to 
Refs. 0, HU for details. 

Finally, whilst the above results and analysis were for the normal hierarchy of neutrino masses, we did also consider 
an inverted hierarchy, where from an analysis of the equations presented in the Appendices we obtain the same features 
and results for neutrino mass runnings (though with different initial values at the Mz scale) . As such, the figures for 
Am^ ol and Am^ tm remain unchanged. Possible changes in the angles and phases arise from the different signs for 
the (nij — mi)/(mj + mi) terms present in each evolution equation, where the 9 12 results remain approximately the 
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FIG. 6: Evolution of 623 in the bulk (left panel) and on the brane (right panel) as a function of the parameter t = ln(^/Mz). 
In this plot we have taken tan/3 = 30. The black line is the MSSM evolution, the red one (small dashes) is for _R _1 ~ 1 TeV, 
the blue (dash-dotted) R^ 1 ~ 4 TeV, and the green (large dashes) R~ 1 ~ 15 TeV. 
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FIG. 7: Evolution of the phase 8 as a function of the parameter t — \n(fi/Mz) with matter fields in the bulk. In this plot we 
have taken tan/? = 30 ( left panel) and tan/3 = 50 (right panel) respectively. The black line is the MSSM evolution, the red 
one (small dashes) is for R^ 1 ~ 1 TeV, the blue (dash-dotted) R^ 1 ~ 4 TeV, and the green (large dashes) i?" 1 ~ 15 TeV. The 
evolution is towards large values at high energies, the divergence signaling the breakdown of the effective theory. 

same (given the relative ordering of masses in these two hierarchies), and the small runnings of #13 and 623 would be 
up rather than down (though as already discussed, these runnings are quite small). 



CONCLUSION 



In the present work we studied the neutrino sector in a minimal supersymmetric model with an extra-dimension 
using the updated value of 6*13 and other mixing parameters. We have presented the evolution equations for the 
physical mixing angles, phases, Am^ ol and Am^ tm , within two distinct scenarios, where a larger tan/3 typically leads 
to more significant rcnormalization group corrections. In the first case, all fields associated with the SM fermions 
are allowed to propagate in the bulk, whereas in the second, the matter fields are restricted to the brane. Neutrino 
masses evolve differently in the two models due to the sign of the (different) dominant contributions in the bulk 
and in the brane cases. For the brane case we find the approximate degenerate neutrino mass spectrum becomes 
more favourable at the ultraviolet cut-off, whilst in the bulk case, the neutrino splitting becomes even more severe 
as the unitarity bounds of the effective theory are rapidly reached. As the evaluation of RGEs may play a crucial 
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FIG. 8: Evolution of the phase 5 as a function of the parameter t = la(fj,/Mz) with matter fields constrained to the brane. In 
this plot we have taken tan j3 = 30 ( left panel) and tan /3 = 50 (right panel) respectively. The black line is the MSSM evolution, 
the red one (small dashes) is for R~ 1 ~ 1 TeV, the blue (dash-dotted) ~ 4 TeV, and the green (large dashes) i?" 1 ~ 15 
Te V. The evolution is towards large values at high energies, the divergence signaling the breakdown of the effective theory. 



role in searching for realistic mixing patterns we also studied the evolution of mixing angles and phases. Contrary 
to the large renormalization effect of #12, the runnings of #13 and #23 were relatively mild. We found a non-zero 
value for #13 during the evolution, which has no appreciable RGE running effects, even when power law evolution 
effects are considered. Therefore it is necessary to introduce new physics effects in order to achieve the tri-bimaximal 
pattern. Here, we also find the maximum CP violation case, 5=5, could be achieved starting from a relatively 
small initial value. In this regard, radiative effects have a very significant impact on neutrino physics. A non-zero 
Jarlskog invariant, which measures the magnitude of leptonic CP violation and is expected to be measured in future 
long baseline neutrino oscillation experiments, could open the door for CP violation in the lepton sector. 
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Appendix A: Conventions for masses and mixing parameters 

Within this appendix we would like to stress our conventions for the mixing angles and phases and briefly discuss 
different scenarios for neutrino masses. The mixing matrix which relates gauge and mass eigenstates is defined to 
diagonalise the neutrino mass matrix in the basis where the charged lepton mass matrix is diagonal. It is usually 
parameterised as follows (26j : 

/ Ci 2 ci3 s 12 c 13 s 13 e- lS \ I e^ 1 \ 

U = -S12C23 - c 12 s 2 3'Si3e 1 ' 5 C12C23 - si2S 23 sx 3 e lS s 23 c 13 e l<fa , (Al) 

\ S12S23 - Ci2C23Sl 3 e il5 -C12S23 - S 12 C 2 3S 13 e tS C23C13 /\ 1 / 

with dj — cos Oij and Sy = sin6ij (ij — 12, 13, 23). We follow the conventions of Ref. [l5[ to extract mixing parameters 
from the PMNS matrix. 

Experimental information on neutrino mixing parameters and masses is obtained mainly from oscillation exper- 
iments. In general Ara^ im is assigned to a mass squared difference between v 3 and V2, whereas Am^ ol to a mass 
squared difference between z^ 2 and v\ . The current observational values are summarised in Table [I] Data indicates 
that Am^ <C Am^ (m , but the masses themselves are not determined. In this work we have adopted the masses of the 
neutrinos at the Mz scale as mi = 0.1 eV, ni2 = 0.100379 eV, and 7713 = 0.11183 eV, as the normal hierarchy (whilst 
any reference to an inverted hierarchy would refer to 7713 =0.1 eV, with m 3 < mi < m,2 and satisfying the above 
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bounds). For the purpose of illustration, we choose values for the angles and phases as the Mz scale as: 612 = 34°, 
13 = 8.83°, 6> 23 = 46°, 5 = 30°, fa = 80° and fa = 70°. 



Parameter 


Value (90% CL) 


sin 2 (26»i 2 ) 


0.861(t°;°i) 


sin 2 (26*23 ) 


> 0.92 


sin 2 (26»i 3 ) 


0.092 ±0.017 




(7.59 ±0.21) x 10 _5 eV 2 


Am 2 tm 


(2.43 ±0.13) x 10~ 3 eV 2 



TABLE I: Present limits on neutrino masses and mixing parameters used in the text. Data is taken from Ref. 0/ for sin 2 (2^13), 
and from Ref. 



Appendix B: Yukawa couplings 

In the followin g w e write for completeness all the evolution equations, where we shall use a notation similar to the 
ones of Refs. [l5lll6j. Note that the beta functions contain terms quadratic in the cut-off, where this part dominates 
the evolution of the Yukawa couplings and of k. The top Yukawa coupling becomes non-perturbative before the gauge 
couplings thus limiting the range of validity of the effective theory. We shall first write down the result of the MSSM 
and then generalize it to the inclusion of the effects arising from the extra dimensional degrees of freedom. 

As in our previous works, the initial values we shall adopt at the Mz scale are: for the gauge couplings we have 
ai(M z ) = 0.01696, a 2 (M z ) = 0.03377, and a 3 (M z ) = 0.1184, and for the fermion masses m u (M z ) = 1.27 MeV, 
m c {M z ) = 0.619 GeV, m t (M z ) = 171.7 GeV, m d (M z ) = 2.90 MeV, m s (M z ) = 55 MeV, m b (M z ) = 2.89 GeV, 
m e {M z ) = 0.48657 MeV, m^{M z ) = 102.718 MeV, and m T (M z ) = 1746.24 MeV [J|28|. 

1. MSSM 

The evolution equations for the MSSM are a limiting case of the 5D models we shall consider in the following. 
In any case, when < t < ln( M * R ) (that the energy we consider for the evolution is Mz < fi < 1/R) the Yukawa 
evolution equations are dictated by the usual MSSM: 

2 dYd _ W ,m„, v t V \ , TJvfvl 1 Qvtv. 1 vfvi v. I J_„2 , o 2 , 16 2 



16^ 2 ^ = Y d (3Tr(YjY d ) + Tr(Y r }Y e ) + 3YJ Y d ± Y^Y U ) - Y d ^-g\ ± 3.g 2 2 ± -g 2 3 
167^ = Y u (3Tr(Y^Y u ) + 3Y„+Y U ± YjY d ) Y u ^gj + ig\ + fg^ , (Bl) 
16tt 2 ^ = Y e (3Tr(Y^Y d ) + Tr(Y^Y e ) ± 3Y e tY e ) - Y e {^g\ ± 3g 



These equations are modified when we enter the energy regime where the effects of the extra dimensions sets in. The 
modifications depend on the particles non-decoupled at a certain energy and on the structure of the model. We shall 
consider two cases, one in which all particles can propagate in the extra dimensions (bulk scenario) and the other in 
which SM particles are constrained to the brane (brane scenario). 

2. The Bulk Scenario 

The RGEs for the Yukawa couplings in the 5D MSSM, for all three generations propagating in the bulk, can be 
expressed as: 

16^ 2 ^ = Y d (3Tr(Y^Y d ) + Tr{Y}Y e ) + 3Y^Y d ± Y^Y u )nS(t) 2 - Y d (^g\ + Zg\ ± fgA S(t) , 
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16tt 2 ^ = Y u (3Tr(Y^Y u ) + 3Y^Y U + YjY d )nS(t) 2 Y u (^g 2 + 3g 2 + fg^j S(t) , (B2) 

16tt 2 ^ = Y e (3Tr(YjY d ) + Tr(Y^Y e ) + 3Y^Y e )-KS{tf Y e f^-g 2 + 3.g 2 ) S{t) 

That is, when the energy scale /z > 1/R or when the energy scale parameter t > ln( M * fl ) (where we have set Mz as 
the renormalization point, and use S(t) — e t MzR)- We can convert these equations to the following form: 

16tt 2 ^ = Y d {T d irS 2 -G d + (3Y^Y d + YjY u )irS 2 } , 

16^ 2 ^ = Y u {t u itS 2 -G u + (3YjY u + YjY d )7TS 2 } , (B3) 
16^ 2 ^ = Y e {T^S 2 -G e + CdY}Y e )nS 2 } , 



where 



T d = 3 Tr(YjY d ) + Tr(Y*Y e ) , 
T u = 3 Tr(Y^Y u ) , 

G u = (^g 2 + 3 9 2 + ^-g 2 )s(t), (B4) 



.r 1 JZ 3 

G e = Qs? + 3s§W). 



3. The Brane Scenario 



The results for the beta function in the case where all matter supcrfields are constrained to live on the 4D brane, 
the quadratic evolution due to the sum over two Kaluza-Klein towers will be milder. The evolution equations for the 
Yukawa couplings are given by: 

16?r2 ^ = Y ^ Tr ^) + T ^Ye) + (QYjY d + 2Y^Y u )S(t)) Y d (^g 2 + \g 2 + |.g 2 ) S(t) , 

ie^ = y u (3Tr(yjy„) + (6yjy„ + iY}Y d )S{t)) y u (^ g { + \ g \ + | ff 2 ) s(t) , (b 5 ) 
16^ = Y e {3Tr{Y}Y d ) + Tr{Y^Y e ) + &Y^Y e s{t))-Y e (^ g 2 + ^ 

Which we can convert to the following: 

16tt 2 ^ = Y d [T d - G d + (6YjY d + 2YJY u )s} , 



16n 2dYu 



Y u [T u -G u + (6Y^Y u + 2Y}Y d )s) , (B6) 
AY 

16^ 2 — - - Y e {T e -G e + (6YX)S} , 



dt 
where 

T d = 3Tr(YjY d )+Tr(YjY e ) , 

/19 2 9 2 32 2 \ 
G d = [ Yo 9i + - 2 9 2 2 + Y 9l)S(t) 

T u = 3 TrY^Y u , 
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Gu = (frf + \gl + f ff3 2 ) S(t) , (B7) 
T e = 3Tr(YjY d )+Tr(Y^Y e ) , 
Ge = (|.9? + \ 9 l) S(t) . 

Appendix C: Mass and mixing angles evolution 

As such, for the four dimensional MSSM, where C = 1 and a = 6Tr(YjY u )- %g\-<og\ = 6(y?+y 2 +y 2 )- §# 2 -6gl, 
we have: 

16tt 2 ^- = mi {a + Cy 2 T (2 sin 2 (6»i 2 ) sin 2 (6> 23 ) - sin(6>i 3 ) sin(26>i 2 ) sin(26> 23 ) cos((5) 

+2 sin 2 (# 13 ) cos 2 (0 12 ) cos 2 (0 23 )) } , 

167r2 ^T = m2 ^ a + Cy * ( 2 cos2 ( 6 ' 12 ) sin2 ( 6 ™) + s H 9 i3) sin(26>i 2 ) sin(20 23 ) cos(J) 



2sin 2 (0 13 )sin 2 (0 12 )cos 2 (0 23 ))} 



16 ^ 2 ^? = TO 3{a + Cy 2 (2cos 2 (0 13 )cos 2 (0 23 ))} , (CI) 
and for the mixing angles: 

2 d9u 1 „ 2 / m i + m 2 ,<Pi — <h\fl 1 4>i — 4>2 



2 m 1 +m 2 ,<P\-<p2,(t ,<Pi -<P2, 2ia \ ■ too \ 

r n o Cy A ZT C0S ( o Ho cos ( o ) cos #13 sin(26>i2 

at 2 \ vfi\ — ui2 2 \2 2 

+ {\ cos( — 2 ^ X~ 3 + cos ( 26l i3)) cos(26> 23 ) sin(2(9 12 ) - [cos <5cos( ^ 1 ~ ^ 2 ) cos(26> 12 ) 
+ sin 5 sin( ^ ^ 2 )] sin(0i 3 ) sin(26> 23 )} 



+ + m 3 
mi — m 3 



— — - cos(<5 - — -) sin(6»i 2 ) sin(0 13 ) ( - cos(— -) sin(0 12 ) sin(26> 23 ) 
— ?n 3 z \ z 

+ 2 cos(<5 - y ) cos(0 12 ) sin(0 13 ) cos 2 (0 23 )) 

' j3 cos(<5- y)cos(0i 2 )tan(0i 3 )(cos(y)cos(0i 2 ) cos(#i 3 ) sin(20 2 3) 



+ 



m 2 — m 3 



+ co S (S - ^ ) sin(0 12 ) sin(2# 13 ) cos 2 (# 23 ))) , (C2) 
167r2 ^T = ^y 2 f^i±^cos(«5-^)cos(0 12 )cos(0 13 )(-cos(^)sin(0 12 )sin(20 23 ) 

at 2 \TTli — 77l 3 2 \ 2 

+ 2 cos(<5 - y ) cos(0 12 ) sm(6 13 ) cos 2 (0 23 )) 

n3 cos(<5 - y) sin(6» 12 )^ cos( y) cos(^i 2 ) cos(6» 13 ) sin(20 23 ) 



m 2 - m 3 



+ cos(S - ^) sin(0 12 ) sin(2# 13 ) cos 2 (# 23 ))^ , (C3) 
16tt — = -Cy T I ^— — cos(y)sin(0i 2 )^cos(y)sin(0i 2 )sin(20 2 3) 
-2 cos(<5 - y ) cos(0i 2 ) sin(0i 3 ) cos 2 (# 23 )) 

+ m2 + m3 cos(%) cos(6»i 2 ) sec(6>i 3 ) f cos(^) cos(#i 2 ) cos(6»i 3 ) sin(20 23 ) 
m 2 — m 3 2 V 2 

+ cos(<S - ^) sin(0 12 ) sin(2# 13 ) cos 2 (# 23 ))) , (C4) 
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As we have already mentioned, the transition to the bulk case will be done by making the replacement of C = 
C(fi) = TTfi 2 R 2 = 7r5(t) 2 and a = QirS(t) 2 Tr(Y^Y u ) — (^g 2 + 6g 2 )S(t). Similarly, we will also have the same equations 
in the brane case, and with C = C(p) = 2^R = 2S(t) and a = 6Tr(YjY u ) - (ffir 2 + 9g%)S(t). 



Appendix D: Dirac and Majorana phases evolution 

As already discussed, in the MSSM C = 1, in the 5D brane case C — 2\iR = 2S(t), and in the bulk case 
C = C{ji) = nfi 2 R 2 = irS(t) 2 . In which case, the equations for the evolution of the phases can be expressed as: 



16tt 



2 aq>\ 1 ^ 2 ( 1 mi + m 2 . ,<fii — (j>2 
~dT ~ 2 



^(2m7^ sin( ^ )cot M 



2 cos(^— — ) cos 2 (6>i 3 ) sin(20i 2 ) 



| C0S (^L_i 2 .)(_ 3 + cos(20 13 )) cos(20 23 ) sin(20 12 ) 



cos 5 cos (- 



'-) cos(2#i 2 ) + sin<5sin(- 



sin(0 13 )sin(2023)}) 



m 1 +m 3 i 
mi — m 3 



2 cos(6>i 3 ) ( - cot(<9 23 ) sin( y) sin(<9 12 ) sec(6>i 3 ) + 2 cos(<9 12 ) sin(<5 - ^-) tan(6>i 3 ) 
+ sec(0i 3 ) tan(0 23 ) sin( y ) sin(0i 2 )) ( - 2 cos(<5 - ^-) cos(6>i 2 ) sin(0i 3 ) cos 2 (6> 23 ) + cos( y ) sin(0i 2 ) sin(2# 23 ) 

+ 21-21 ( sin{5 ~ ^ sin(012) tan(013) + cos ^ 12 H cot ^ 23 ) sec ^ 13 ) sin (y) 

— cos(y) cot(#i 2 ) sin<5tan(^i 3 ) + sin(^) cos($i 2 ) cos S tan(#i 3 ) — sec(#i 3 ) tan(# 23 ) sm (^)}) 

x (-2cos(^)cos(0 12 )cos(0 13 )sin(20 23 )-2cos( ( 5-^)sin(0 1 2)sin(20 13 )cos 2 (02 3 ))j , (Dl) 



16tt 



2 rf^2 ^ /-y„.2 n m l+™2 ... 



= -Cy 

2 T V 2 mi — m 2 



( ^L_ia } tan( 12) ( _ 2 cos( ^L_i 2 .) cos 2 (0i 3 ) sin(20i 2 ) 



{cos(- 



-)(-3 + cos(26»i 3 )) cos(2(9 23 ) sin(26>i 2 ) - 4[cos<5cos( — 



-)cos(20i 2 ) 



+ sin^sin( 01 02 )]sin(g 13 )sin(2g 23 )}) 

- Wl + m3 2cos(6>i 3 ) f cot(0 23 ) sin(^-) sin(0i 2 ) scc(#i 3 ) - cos(6»i 2 ) sin(5 - ^-) tan(#i 3 ) 
mi — m 3 V 2 2 

+ sin(6*i 2 ){sin(5cos(^-) tan(#i 2 ) tan(#i 3 ) — sin(^-) cos<5 tan($i 2 ) tan(#i 3 ) — sin(^) scc($i 3 ) tan(0 23 )}^ 
x ( - 2 cos(<5 - ^) cos(0i 2 ) sm(6 13 ) cos 2 (fl 23 ) + cos(^) sm(6 12 ) sin(2# 23 )) 

+ m2 + TO ' 3 f sin(^) cos(^i 2 ) cos(2# 23 ) csc(6> 23 ) sec(6»i 3 ) scc(6> 23 ) + 2 sin(c5 - ^) sin(0i 2 ) tan(#i 3 )) 
m 2 — m 3 \ 2 2 / 



x ( - 2 cos(^) cos(0 12 ) cos(0i 3 ) sin(2# 23 ) - 2 cos(<5 - ^) sin(0 12 ) sin(2# 13 ) cos 2 (0 23 )) 



(D2) 



1 a 2 dS 1 ^ 2 ( m l + m 2 . >1 ~ 

107T — = -6V sm( — 



; -)csc(20i 2 )( - 2cos(- 



; -)cos 2 (6li 3 )sin(2i9i2) 



{ 



cos 



— )(-3 + cos(26»i 3 )) cos(26» 23 ) sin(20 12 ) - 4[cos5 cos(- 1 



2 " v "' v " v 2 

+ sin(5sin( ^ 1 - ^ 2 )] sin(gi 3 ) sin(2g 23 )|^) 

- g^t^) ( cos ( y){ 2 - 6 cos(2#i 2 ) + cos(20i 2 - 20i 3 ) - 6 cos(2#i 3 ) 



-)cos(2fl 12 ) 
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cos(26»i2 + 20i 3 )} csc(0i 3 ) sec(6» 12 ) sin ,5 + 4 cos(6> 12 )(-3 + cos(26>i 3 )) csc(6>i 3 ) sm(S - — ) 
sin(y)| - cos(5)[2 - 6cos(26» 12 ) + cos(26> 12 - 26> 13 ) - 6cos(26> 13 ) + cos(26> 12 + 20 13 )] csc(0 13 ) sec(0 ; 
16 sin(0 12 ) cot(0 23 ) - 16 sin(0 12 ) tan(0 23 )}) ( - 2 cos(<S - y ) cos(0 12 ) sin(0 13 ) cos 2 (0 23 ) 
cos (y ) sin(6»i 2 ) sin(20 23 )) 



+ cos(y)sin(0i 2 )sin(20 23 ) 

" 4(m 2 2 -m 3 3 ) ( [ 1T + I COS ( 2012 ) " i cos ( 20 i2 - 2 M ~ 5 cos(20 13 ) - ~ cos(20 12 + 20 13 )] 
x csc(#i2) csc(^i3) sec(^i3) sin((5 — — ) + cos(#i2) sec(#i3)[— cos( — )(— 3 + cos(2#i3)) cot(#i2) csc(0 
;(0 23 )+sm(y^ 



2 

- 4 sin(y ) cot(6» 23 ) + sin(y ) cos 8{-Z + cos(20i 3 )) cot(6>i 2 ) csc(0i 3 ) + 4 sin(y ) tan(0 23 )] 
x ( - 2 C os(y ) cos(0i 2 ) cos(0 13 ) sin(20 23 ) - 2 cos(<5 - ^) sin(0 12 ) sin(20 13 ) cos 2 (0 23 ))^) . 



'i 3 ) sin 5 



(D3) 
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